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ABSTRACT 



We report a technique to calculate the impact of distinct physical processes inducing 
non-Gaussianity on the cosmological density field. A natural decomposition of the cos- 
mic genus statistic into an orthogonal polynomial sequence allows complete expression 
of the scale-dependent evolution of the topology of large-scale structure, in which effects 
including galaxy bias, nonlinear gravitational evolution and primordial non-Gaussianity 
may be delineated. The relationship of this decomposition to previous methods for 
analysing the genus statistic is briefly considered and the following applications are 
made: i) the expression of certain systematics affecting topological measurements; ii) 
the quantification of broad deformations from Gaussianity that appear in the genus 
statistic as measured in the Horizon Run simulation; iii) the study of the evolution 
of the genus curve for simulations with primordial non-Gaussianity. These advances 
improve the treatment of flux-limited galaxy catalogues for use with this measurement 
and further the use of the genus statistic as a tool for exploring non-Gaussianity. 



Subject headings: cosmology: theory — dark matter — methods: statistical 



1. Introduction 

The study of the topology of large-scale structure promises much from the current generation 
of galaxy redshift surveys. On the largest scales, the relationship between the development of struc- 
ture and primordial non-Gaussianity has recently deepened; in the regime of nonlinear structure 
formation, the effects of gravitational evolution and the complicated phenomenon of bias continue 
to stimulate interest in the flelds of cosmology and galactic astrophysics, where, however, the use 
of the (geometric) correlation function still predominates. The reason for this seems to be that the 
topology of large-scale structur e is difficult to interpret when the underlying distribution departs 



Ghoi et al.ll2ninl ). 



from a Gaussian rand om field (jPark et al.ll2005l : I James et al.ll200?i : iGott et al.ll2008l : iJames et al 
20091: 
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Of the four Minkowski functionals t hat ch aracterize the morphology of the cosmological density 
field, the genus statistic of iGott et alj (|l986l ) has been most frequently employed. Given a choice 
of critical density value, the genus is a measure of connection or isolation of structures enclosed by 
the surface of that constant critical density. Different choices of density threshold excise filaments, 
cluster and voids, whose topological physics is then expressed in compact form. The curve of genus 
as a function of density has a known analytical form when the underlying density field is of Gaussian 
random form — the task is to understand the manifestation of non-Gaussian physics in the topology 
of cosmic structure. 



1.1. Regarding density thresholds 



To help with the development of these ideas, several prevailing definitions of density thresh- 
olds are reconciled. The parameterisation of density used in genus measurements is idiosyncratic, 
but well-motivated; it is the simplest mapping independent of the underlying density distribution 
function. There are two (less preferable) alternatives: for density fields 6 = p/p — 1 whose dis- 
tribution functions f{6) depart from Gaussian form — as the cosmological density field must over 
at least some scales — the parameter u = 5 /ap is a misleading mapping, imposing as it does an 
artificial sy mmetry about r nean d ensity. The often-mentioned lognormal model for the distribution 
function of IColes Joned (Il99ll ) improves this insofar as it corresponds to the distribution func- 
tion over a broader range of scales: using the parameter 6 = ln(l -|- S), the mapping is given by 
= 5/ag. Except in highly evolved fields, this mapping is nearly equivalent to the most general 
parameterisation: 



where fci^) = Erf 



fi5')d6' 



(1) 



and it is this variable that represents density in topological statistics. This choice of mapping insists 
that regions should be judged as comparably over- and under-dense when they occupy the same 
volume, rather than on the basis of their density. 



2. A Hermite function transform 

The genus of a two-dimensional surface is a measure of its connectedness: the number of 
holes through the surface less the number of isolated regions (plus 1, by convention); a donut 
has genus 1, a sphere and two spheres —1. The genus curve of a Gaussian random field, and 
that of the primordial U niverse in the simplest models of inflation, is understood to have the 



form (jBardeen et al.lll983l ) 



5GRf(i^) oc exp ( -y ) (1 - , (2) 
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with the constant of proportionahty a functional of the power spectrum of the field; in any event, 
it is the shape of the curve that is to be scrutinized. 

The curve of genus number as a function of density provides information even in this extremely 
restricted case: structures separating the highest-density half of the Universe from the lowest- 
density half are many times connected; excising extreme densities yields many disjoint regions. 
The balance between these two regimes and the relative abundance of isolated high and low density 
regions are a scale-dependent probe of many physical processes. Consequently it is desirable to 
study the surface of genus number as a function of density and scale. In practice, different scales 
are analysed separately before being compared. 

To address the manner in which distinct physical processes alter the topology of large-scale 
structure, it is proposed that the genus surface (i.e., the genus curve as a function of scale, char- 
acterized by Gaussian smoothing length A) be decomposed in an orthogonal basis of Hermite 
functions: 

oo 

9{^;X) = ^anWi^ni^) (3) 

n=0 

/oo 
g{u]X)ipn{y)dv. (4) 
-oo 

The Hermite functions il^niy) are themselves weighted analogues of the Hermite polynomial^ 



n 

2 



e--/2 (5) 



= ^ e--''^Hn{y) (6) 
V n! v27r 

which will satisfy the orthogonality relation 

lljm{y)i^n{y)dv = 5mn (7) 



required for the decomposition @ to exist. The evolution of the Hermite modes encode, as a 
function of scale, the imprint of the physical processes that have modified the field from the Gaussian 
random form. The sequence of coefficients a„ is referred to as the Hermite spectrum of the genus 
curve. The normalised sequence of coefficients, o„ = |an|/ Xlm 1'^"^! employed throughout 

this work as it provides a satisfying interpretation of the fractional contribution of each Hermite 
mode as a probability distribution across the modes, with different physical processes mapping 
onto different distributions. In particular, the quantity 02 is a completely general measure of the 



^There are two definitions for the Hermite polynomials, conventionally referred to as the 'physicist' and 'probabilist' 
forms. The probabilists' definition is used in this communication, as the genus curve for a Gaussian random field is 
more naturally expressed this way. 



-4- 



Gaussianity in the field: for a Gaussian random field, 02 = 1; for all non-Gaussian fields, 02 < 1, 
where the particulars of the non-Gaussian process are encoded in the other coefficients. 

Much of the arithmetic is simplified by the rescaling g{u) — >■ e^^/^g{u), which ensures that 
any Gaussian random field has a trivial decomposition: in this scheme, the base mode of the 
transform is n = 2, i.e., gcKeW) oc ^2(^^)1 with other low modes appearing as the field is perturbed 
from a Gaussian state. Odd-numbered modes will, in particular, introduce asymmetric features 
corresponding to an overabundance of either isolated clusters or voids. In practice, this means 
that a measured genus curve should be multiplied by a factor of e^^/'^ before the decomposition, a 
convention adopted throughout this work. 

Quantifying the distortions to the genus curve due to a physical process inducing non-Gaussianity 
has long been understood to deserve attention. Two approaches that have led to progress are: i) 
the analytical method invoking second-order perturbation theory to calculate the impact of small 
departures from Gaussianity; and ii) the use of a phenomenological set of gcnus-rclatcd statistics 
that quantify differences between a measured genus curve and that for a Gaussian random field. In 
the following sections, the formalism described above is applied to the examination of a range of 
physical and unphysical (i.e., systematic effects) phenomena that create non-Gaussian signatures 
in the cosmological density field. 

3. Application I: Formulae for treatment of density field systematics 

In this section, the formalism is applied to the task of correcting observational measurements of 
the genus statistic for systematic selection effects. These different effects distort the genus curve in a 
non-trivial way, to which correction formulae arc derived as a function of the relevant reconstruction 
parameter. These correction formulae make use of the Hermite decomposition technique for genus 
analysis. 

To examine these systematics, use is made of simulated Gaussian random fields of power 
spectrum P{k) oc k^, though this choice of function is arbitrary up to the amplitude of the resulting 
genus curve. The fields are constructed in a spherical polar Fourier space at each wavenumber k 
by assigning Fourier amplitudes from a Rayleigh distribution of parameter P{k), phases from 
a distribution uniform in the interval [0, 27r) and applying the conjugate (Hermitian) symmetry 
required to ensure a real- value configuration-space density field. 

3.1. Sampling systematic 

To return a fair reconstruction of the density field from a density-dependent sampling (i.e., 
nocp^6~l)ofA'' points, it is necessary to smooth the field with a filter exceeding a characteristic 
length scale of the point distribution. When this condition is not satisfied, the density field is not 



- 5 - 



adequately reconstructed at that scale, and the genus measurement of the smoothed field will not be 
accurate. There is excellent motivation to understand and correct this effect, as doing so improves 
the utility of flux-limited samples for the purpose of the genus measurement, offering order-of- 
magnitude improvements in signal relative to a volume-limited-sample. This subsection provides 
a parametric form for this correction and demonstrates its application to the case of a Gaussian 
random field with power-law power spectrum. 

To examine this systematic, the same Gaussian random field realisation is Poisson sampled 
Ur = 100 times. The variance associated with the correction to the statistic is not explored in this 
work. Each sample is smoothed through a Gaussian window of scale parameter A and the genus 
curve for the recovered field is measured. The quantity used here to parameterise the distortion 
present in the genus curve is the ratio of a characteristic length scale in the point distribution to 
the smoothing length: 

T = j; (8) 

the quantity i can admit a variety of definitions, the choice of which is only important in so far as it 
be treated consistently. Among these are the radius of the average spherical volume per point; the 
side length of the average Euclidean volume per point; the empirical average separation between 
every pair of points; and the average nearest-neighbour separation of a Poisson distribution with the 
same sampling density; indeed, it seems that the definitions return length scales that are equivalent 
up to a constant factor. For the purpose of this work, the first of these 



(—) 

\ Attn J 



1/3 

(9) 



is chosen. 



It is hypothesised that the systematic is a function of T, so that an additive correction factor 
to the raw genus curve can be defined 



tmax 



f{u;T)= }_^an{T)M^y, (10) 

n=l 

for a given input power spectrum, the task is to empirically derive the correction coefficients a„(T), 
which at best will require tabulation only to very small rimax- This proves to be the case for power 
spectra that have been tried: the correction function is accurately described once Hermite functions 
up to fourth order are included, with no significant gains beyond this. Furthermore, above values 
of T ~ 1, the correction is in the first Hermite mode at the 95% level, so that a formula can be 
given in terms of this quantity only. A figure showing the coefficient amplitude oi as a function 
of T is given in Figure [TJ Extension to arbitrary T is achieved by interpolation or the use of the 
function: 

ai(T) = /3il + /3o; ^^^Bo^; (11) 
T 02(1) T 
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T = Xjt 

Fig. 1. — Sampling systematic correction factor for a point-sampled Gaussian random field: as 
the smoothing length falls below the mean separation between points, the asymmetry in the genus 
statistic is accurately corrected by the first-order Hermite function with amplitude ai(T) (top 
panel); the magnitude of this correction relative to the amplitude of the base function, the second- 
order Hermite function, is shown as the fractional correction (bottom panel). This example shown 
here is for a field of power-law power spectrum with index n = 2. 
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in the example case presented here, the maximum-hkehhood estimates of the parameters are /3i = 
—0.0548 and /3o = 0.0166, and Bq = 0.13. It is argued here that this formula is robust with respect 
to the choice of power spectrum, absolute sampling density and the resolution of the density field. 
Considering the form of the genus curve presented in Equation (2), it is well-known that the 
power spectrum of the Gaussian random field enters only in the constant of proportionality. For 
small departures from Gaussianity, such as those studied by Matsubara (2003), the genus curve 
is modified by low-order Hermite polynomials, as is the case here, with the coefficients of these 
polynomials depending on the power spectrum via the skewness parameters. Crucially, because the 
Hermite mode coefficients a„ are calculated only with respect to the shape of the curve and not to 
its amplitude, the leading-order infiuence of the power spectrum is removed. 

The result is that the parameterised form of the correction presented in Equation (|lip is 
applicable to a wide variety of power spectra, though the parameters themselves will need to be 
tuned in each case. It is foreseeable that an analytic form of these parameters could be derived in 
the same fashion as the second-order perturbation theory analysis of non-Gaussianity arising from 
gravitational evolution, though this is beyond the scope of the current exposition; an attempt to 
demonstrate both the generality of this form and the derive a connection between the parameters 
and the underlying power spectrum would be a worthy future advance. 

It is clear, however, that this correction is not useful at arbitrarily small T. Indeed, while the 
prevailing rule of thumb that T ~ 1 is sufficient to ensure a reasonable reconstruction seems to 
be borne out, there is nevertheless some correction required even then. While it is possible to go 
somewhat below this, the corrections quickly begin failing to converge; treatment when T < 0.75 
seems ill-advised — better to deem the field insufficiently reconstructed. Above this limit, however, 
the correction formula holds to better than 1%. Further work will clarify the degree to which this 
correction holds for non-Gaussian {e.g., nonlinear) fields, exploring whether the small departures 
in Figure [1] are a breakdown of this scheme. 



An extension of the sampling systematic correction is its application to flux-limited galaxy 
samples in which the survey selection function is well understood. When the selection function is 
known, it can be used to upweight those locations in the density fleld that have been relatively 
under-observed (for whatever reason), prior to smoothing. Nevertheless, these regions will still 
have lower sampling rates, which will induce a sampling systematic that varies with location. To 
take the speciflc example of redshift, at high redshifts the average pair separation i will often be 
relatively higher even with the selection upweighting taken into account. Because the genus is a 
locally additive quantity (after the reconstruction has taken place), the formula ([IT]) can be used 
to correct the resulting measurement: 



3.2. Selection function systematics 




(12) 
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The improvement in signal-to-noise that results from using flux-limited samples depends on the 
specific selection properties of the survey, as the quality of the reconstruction of the density field is 
proportional to the number of galaxies used in its reconstruction; a volume-limited sample typically 
contains at most 10% of the total sample, which propagates directly into the resulting genus signal 
(see, e.g., James et al. 2009). However, the type of sample that should be constructed will ultimately 
depend on the scientific application: a flux-limited sample will, in general, contain greater variation 
of the galaxy mass function with redshift, which will impede the interpretation of measurements. 
The purpose of the study of the sampling systematic is to ensure that, when the selection and 
galaxy sample properties are well understood, Gaussian smoothing reconstruction can be used to 
study scales close to the mean inter-object separation. 



3.3. Edge effect systematic 



The edge systematic has been treated in some depth bv iMelott Sz DominikI (|l993l ). who pi- 
oneered a windowing correction that modifies the reconstructed field before genus measurement, 
rather than correcting the measurement after it is made. The motivation for the existing method is 
that power will be pushed through these edges during the reconstruction, so that boundary regions 
should be expected to be systematically under-dense. This method works by creating a field of 
constant density value with the same geometry as the data, performing the same reconstruction 
procedure on it and ratio the smoothed data field by the smoothed constant field. This remains the 
best correction method for this effect, performing well for the current generation of surveys. These 
authors have noted that even when their correction method is applied, the measured genus curve 
should be scaled by a factor corresponding to the ratio of the total initial volume to the volume 
not pushed outside during the smoothing process, which is carried through to the computation 
of Hermite coefficients. As this affects all coefficients equally, their fractional contributions are 
unaltered. 



3.4. Discretisation systematic 

Though the detailed specification of this systematic effect is unclear, it is known that a min- 
imum smoothing length exceeding the scale of individual pixels must be applied to any pixellated 
density field. When this scale of smoothing is not met, the genus measurement is distorted sym- 



metric ally about median density in the manner shown in Figure [21 In work by [Hamilton et al 



(|l986l ). the effect of finite (octahedral) pixel size was developed analytically, though without ex- 
plicit invocation of a smoothing scale (c/. Equation (21) of that work). This computation shows 
that the expected expression (Equation (2) in the present work) is recovered in the limit of van- 
ishing pixel size r and that residual terms 0{r'^) are well-described by q uantities formulated from 



the higher spectral moments of the power spectrum. iMelott et al.l (|l989l ) have developed this idea 



in the context of the two-dimensional genus statistic, but with application to cubical pixels of the 
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Fig. 2. — Demonstration of the discretisation systematic effect, in which the genus curve of a 
pixellated Gaussian random density field (soHd) is symmetrically distorted from its theoretical form 
(dashed), where each curve has been normalised so that the absolute sum of its Hermite coefficients 
is 1. The top panel shows the difference between the two curves; the right-hand inset shows the 
Hermite spectrum recovered from the decomposition of the measured genus curve, demonstrating 
that the systematic effect is not corrected by the treatment of only low-numbered modes; but, as 
in the left-hand inset, the formula for its correction obeys the relatively simple form of equation 
()13p . shown here with the data (points) and models (solid) for A — )• (no smoothing), A = 0.5 and 
1.0, tending to the limit of a delta function at n = 2, corresponding to the Gaussian random field. 
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kind most commonly employed in genus calculations. This impressive computation, to sixth order 
in the pixel size, shows clearly how polynomials in such as the Hermite polynomials, carried to 
progressively higher orders with coefficients expressed in terms of spectral moments, can account 
for the behaviour observed in the large pixel-size regime of genus computations. 

We present here a treatment that complements the formalisms that have been developed 
previously, both by explicitly mixing consideration of the smoothing scale and by presenting a 
form of that correction that carries to all orders. For this systematic, the bias induced in the 
genus curve is not a function of a single Hermite mode, but of the full Hermite spectrum of the 
measurement. This complicated form (see right-hand inset in figure) suggests that treatment for 
this systematic can be achieved by the inclusion of corrections to Hermite coefficients at all orders. 
The proposed formula for this correction is 

j 0, n = l,3,5,... 

= [ (-ir/^e--" n = 2,4,6,... ' 7- 7^ + 70, (13) 

where the maximum-likelihood parameter values are 71 = 3.75 and 70 = 0.325 for the field presented 
in this example. The left-hand inset in Figure [2] demonstrates the variation of this systematic for 
several smoothing lengths (A = 0, 0.5 and 1.0 pixels, tending toward the delta function 5{n — 2) 
in the limit A — )• 00. Even when A exceeds the pixel size, some effect remains. In practice, this 
systematic is not often an issue because the smoothing scale will naturally exceed the pixellation 
length by a sufficient factor in order to ensure that the sampling systematic is addressed, but further 
work can develop a treatment for cases where this is not achievable. 

A promising direction in which to extend this result is to consider expressing the coefficients of 
the correction in terms of the spectral moments of P{k), as has been done in earlier work. While not 
attempted here, some guidance can be found by comparing the terms of the expansions presented 
in those previous results to the forms of Hermite polynomials, matching coefficients as polynomials 
of spectral moments. 



4. Application II: Astrophysical and cosmological phenomena 



4.1. Nonlinear gravitational evolution 

4.1.1. Perturbation theoretic treatment 

An infl uential forrnalism for the study of the topology a t weakly nonlinea r scales has been 
provided by iMatsubaral (119941 ) and developed energetically by iMatsubaral (120031 ). In this pertur- 
bative scheme, carried out in powers of the variance of the density field, the modes corresponding 
to n = 1,3 and 5 appear with greater relative strength at scales below some scale, A ~ 10 Mpc/h 
at the present epoch. This scheme is naturally re-expressed in terms of the current work, as a 
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prediction for the spectrum of Hermite coefficients, yielding 

5 



(14) 



n=l 



with 



(15) 



5(2)a(A) 

^/65(lV(A) 


vW3)5(°V(A) 

where 5*-'^'^'^^ are the skewness parameters of Matsubaral ( 2003 ) and a is the rms of the evolved 
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density field at spatial scale A. By expressing the non-Gaussian evolution in terms of the relative 
contribution of each mode, the result is independent of the power spectrum of the density field, 
which enters only as the overall normalisation of the genus curve and so scales all modes equally; 
all other coefficients are zero, to this order in the perturbation expansion. 

These two approaches to non-Gaussian effects on cosmic structure topology lead naturally 
to the suggestion that this study be carried out in terms of the Hermite spectrum of the genus 
curve. The ansatz for this programme is that the distinct processes inducing non-Gaussianity in 
the cosmological density field — survey selection effects, bias, redshift space distortions, nonlinear 
gravitational evolution and primordial physics — have differing spectral imprints that will be delin- 
eated. In the remainder of this communication, results applying this decomposition to the study 
of nonlinear evolution and selection effects are presented. 



4-. 1.2. Examination of Horizon Run simulation 
In this se c tion, the formalism of Section 3 is applied to genus measurements in the Horizon Run 



of iKim et al.l (120091 ). a dark matter-only A^-body simulation of some 70 billion particles across a 



v olume of = (6.592 Gpc//i)'^, from which dark matter halos are drawn following the methodology 



of lKim &: ParkI (120061 ). The simulation volume studied here is a periodic cube at redshift zero, with 
halos of all masses down to resolved limit used, though, following convention, each object is assigned 
equal weight, rather than a weight proportional to its mass, during the reconstruction process. 

The topology of structure in this simulation is examined across scales from 5 to 65 Mpc//i. 
At all scales, the spatial resolution of the density field and number of resolution elements in the 
measurement, are held fixed. To achieve this, subvolumes extending from the origin to Lmax = lOOA 
Mpc//i are extracted from the full simulation and embedded within an array of size 
Periodicity is assumed during the density field reconstruction. 

Figure [3] shows the genus curves measured from the Horizon Run subhalo catalogue on scales 
ranging from 5 to 65 Mpc//i. It is surprising to see an asymmetry toward isolated high-density 
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Fig. 3. — Genus curves for the Horizon Run subhalo population on scales of 5 (most black) to 
65 (most blue) Mpc//i, in steps of 5 Mpc//i, normalised to the number of smoothing kernels in 
the total survey volume. Asymmetry in the genus curve, most pronounced at the shortest scales 
where sampling density effects are relevant, nevertheless persists well into the linear regime (see 
Figure HI). The mean inter-halo separation for the Horizon run is 15 Mpc//i; the corrected curve at 
10 Mpc//i may still be deemed useful, while the shortest scale presented in this plot is included as 
a demonstration of the behaviour of asymmetry in the strongly under-sampled limit. 
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structures even at the largest scales of r econstruction, wh ich is well into the linear regime; this 
result confirms the analysis carried out bv lKim et al.l (|2009l ). which is in contrast to earlier studies. 
To explore this measurement further, the Hermite decomposition of the genus curves is presented in 
Figure m in which the contribution from each mode is normalised to give its fractional contribution 
to the total: a„ = lonl/ 1*^™!- '^^^ decomposition expresses concisely the degree of nonlinear 
evolution in the density field as reconstructed from these objects. It is to be expected that part of 
the non-Gaussian signal in the decomposition arises from the use of biased tracers of the density 
field, though this is unlikely to manifest itself in a scale-dependent fashion at the scales studied 
here in the absence of any primordial non-Gaussianity, so that the relative evolution of the odd- 
numbered Hermite modes (Figure H]) is taken to be a description of the modification to topology 
induced by nonlinear gravitational evolution. At the shortest scales, the sampling density falls 
below that required for direct reconstruction of the field, so that the systematic correction derived 
in the previous section is applied to the measurement. Even with this correction applied, it is clear 
that the mode corresponding to a pure Gaussian random field is under-represented relative to the 
expectation from second-order perturbation theory, particularly evident at the largest scales. 



4.2. Redshift-space distortions 



The peculiar motion of galaxies affects the genus measurement only through its impact on 
the reconstructed density field. No longer a systematic to be overcome, this coupling between the 
density and velocity fields yields information about the processes of gravity and structure formation. 
The linear treatment of this coupling, the Kaiser limit, has been shown to generate no distortion 
to the shape of the genus measurement for a Gaussian random field; like the survey boundary 
correction, this affects the amp litude of th e genu s curve, so that the fractional contribution from 
Hermite modes is unchanged. iMatsubaral (jl996l ) has speculated that the short scale distortion 
effects, which act orthogonally to those in the Kaiser limit, w ill also alter just the amplitude of the 
curve, a claim that has been examined in some simulations (jPark et al.ll2005l ). Given the interest 
directed toward redshift-space phenomena in the field at present, it seems plausible that in the 
intermediate limit between these two regimes could host effects that do alter the shape of the genus 
curve, though it is not clear that such scales can be examined in the current generation of redshift 
surveys. 



4.3. Primordial non-Gaussianity 

Broad classes of infiationary models generate significant non-Gaussianity in the early-time cos- 
mological density field. The customary formalism for this effect is through the (Bardeen) curvature 
perturbation in the matter era: 

$ = + /nl(</'' -(</>>'), 
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A (Mpc/h) : Structure scale 



Fig. 4. — Hermite spectrum for the Horizon Run genus curves in Figure [3] as a function of scale: 
fractional contribution of modes o„ from the nonlinear to linear regime (points and solid lines, 
with the mode n = 1 in black), compared to the predictions from second-order perturbation theory 
(dashed). The sampling correction of equation (jlip has been applied to the n = 1 mode; the 
uncorrected mode is shown in grey. 
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where (/> is a purely Gaussian random variable and /nt , is taken to characterise the magnitude of 
non-Gaussianity in the field (jKomatsu &: Spergell 120011 ). The bispectrum of the primordial cur- 
vature perturbation, identically zero in a Gaussian random field, will peak in regions of Fourier 
space for which the closed triplet of wavenumbers {ki,k2,k3} exhibits one of a handful of ratios, 
giving rise to a classification of primordial non-Gaussianity in terms of triangle configurations. 
Of concern here are the 'squeezed' (c/. 'local') {ki ^ A:2 ~ k^), equilateral (/ci ~ ^2 ~ ^3) and 
enfolded {ki ~ A;2 ~ ^3/2) configurations, corresponding respectively to inflations scenarios involv- 
ing mul tiple fields, horizon crossing and a pproximate shift symmetry in very general single-field 
models (jWagner. Verde &: Boubekeuiil20ld ). 



A study of the influence of primordial non-Gaussianit y on the topology of large-scale st r uctur e 
has been undertaken on the scales above 100 Mpc/h by iHikage. Komatsu &: Matsubaral ()2006l ). 
who present formulae for the genus statistic (among others) from perturbation theory. In com- 
paring these analytic results to iV-body simulations, they note that the impact of primordial non- 
Gaussianity should be slight in proportion to the degree o f gravitationally induced non-Gaussianity 
in the field. Nevertheless, the discovery by iDalal et al.l ()2008l ) of a scale-dependent galaxy bias 
induced by local-type primordial non-Gaussianity provides sufficient motivation to attempt an ex- 
ploration of this effect through the formalism presented in this work. 

Here, the impact of primordial non-Gaussianity on the t opology of the late-time density fiel d 
has been evaluated empirically with the simulation data set of lWagner. Verde fc Boubekeurl ()2010l ). 
These data are the output of a sequence of A^-body simulations with generic non-Gaussian con- 
ditions, where the departure from Gaussianity is characterised completely by the bispectrum of 
the initial conditions. Power spectr a, halo mass functions and bias pr opert ies of these sim- 
ulatio n data have been presented by IWagner. Verde fc Boubekeurl (j2010l ) and IWagner &: Verde 
(j201ll ). In this suite, halo catalogues derived from an evolved density with initial non-Gaussianity 
/nl = {0, ±100, ±250, ±500} have been produced, in each of the local, equilateral and enfolded 
configurations. Snapshots have been provided at z ~ {0, 0.33, 0.67, 1, 1.5, 2}; in this work, we con- 
sider halo catalogues from the redshift one density field with a view to optimizing the number 
density of halos and (reciprocally) the degree of nonlinear evolution in the underlying density field. 
Concurrent work by Wales, James &: Lewis (2012) incorporates an examination of the redshift 
dependence of these measurements. 

The halo catalogue output of these simulations is transformed into a counts-in-cells density 
field in the same manner as the data in the previous section. To evaluate the departures from non- 
Gaussianity present in the genus measurement of these data, the magnitude of the base (n = 2) 
mode of the decomposition, relative to that of the Gaussian case, is studied across the range of 
scales from 15 to 65 Mpc//i, in steps of 5 Mpc. The scale at which T = 1 is around 12 Mpc, so all 
scales examined are free from the impact of the sampling systematic. Genus curves are measured 
for each (/nL; configuration) pairing, in each of the eight available simulations, before the Hermite 
decomposition is performed on each measurement. The variance between the measurements of the 
Hermite coefficients for the eight simulations is taken to be an estimator of the variance of the 



- 16 - 



coefficient statistic. 

Figure [5] shows the relationship between /nl and the relative-to-Gaussian strength of the n = 2 
Hermite mode. This is an interesting result that has several possible interpretations. Perhaps 
the most satisfactory is that, as the amplitude of the genus curve is a functional of the (biased) 
power spectrum, even though the underlying dark matter power spectrum is fixed for a given 
redshift, configuration and realisation, the scale dependence of the quantity a2(/NL)/a2 (gauss) is 
a measure of scale-dependent bias. A prediction of this hypothesis is that all other such ratios 
On (/NL)/an (gauss) should show the same effect, however these data do not allow such a test to 
be carried out at sufficient significance; e.g.. Figure [6|) shows the same measurement as Figure [5] 
with the lower significance modes of the genus curve. In this hypothesis, it is significant that 
the local configuration shows this scale dependence. Figure [7] shows how the other configurations 
of the bispectrum do not result in the same behaviour. In some models it is predicted that the 
equilateral configuration should show a 1/k scale-dependent bias, which would be expected to 
manifest here. Even though this is weaker than the l/k"^ dependence of the scale-dependent bias 
in the local configuration, it is somewhat surprising that only the lo cal configuration displays this 



behav iour — however, this seems to be compatible with the results of lWagner. Verde &: Boubekeur 



mm- 



An alternative explanation that is tractable by the means presented in this work is that, if the 
bispectrum manifests in the genus curve in a manner that is not simply related to the amplitude, 
then the results in Figure [8] constrain the mechanism of this action. It might be hoped that, given 
the exact prediction for the genus curve of a Gaussian random field, many distinct types of non- 
Gaussianity could be examined, even phenomenologically. This is a topic that requires a wider set 
of simulation data to examine and that will be explored in the future. 

To return briefiy to the hypothesis of scale-dependent bias manifesting in the amplitude of the 
genus curve: recent work by several authors has given rise to met hods of constraining dark energy 



and modified gravity by means of the genus curve amplitude (IWang et al.l 120101 : IZunckel et al 



2011 



Park &: Kim 



based on the development of a standard ruler from large-scale structure by 
(|2010l ). This class of tests probes the cosmological parameters via the linear regime density field as 
well as the comoving volume element, both of which determine the genus amplitude. Cosmological 
density fields that possess some primordial non-Gaussianity will also preserve their initial (non- 
Gaussian) topology into the linear regime. The impact of primordial non-Gaussianity, as measured 
and interpreted in the present work, can be used as part of the standard ruler technique, by 
exposing an additional scale-dependent term (that is, beyond the usual scale dependence of the 
genus curve amplitude) that does not manifest within the comoving volume element. This improves 
the possibility of distinguishing these phenomena using the genus curve. 
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Fig. 5. — The scale dependence of the 02 coefficient of the genus curve for varying /nl in the 
local configuration relative to the Gaussian (/nl) case, averaged over eight simulation realisations 
(solid), with la error bars shown. The points are displayed in a staggered manner at each scale for 
clarity. The maximum likelihood fit of a one-parameter model of linear scale dependence (dashed) 
indicate a relationship between the rate of scale dependence and /nl- 
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Fig. 6. — Empirical demonstration of the diminshed impact of /nl primordial non-Gaussianity on 
the weaker modes (n = 1 and n = 3) of the redshift zero genus curve, a result that appears to hold 
across the scales available in the simulation set used in this work, with the same labelling scheme 
as Figure [5l The diagonal elements of the covariance matrix, shown as la error bars, have been 
offset for clarity at each scale. 
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5. Discussion 

The challenge posed by this result is to delineate distinct physical processes in terms of their 
impact on particular Hermite modes (or, more generally, sequences of such modes). This article has 
introduced a Hermite decomposition of the genus curve for cosmic structure, Eq. ^ that serves as a 
pivot between genus measurements, physical modelling such as that from second-order perturbation 
theory and the systematic effects that arise in the measurement of the genus statistic in galaxy 
redshift surveys. In contrast to existing genus-related statistics, because the set of coefficients 
ttn is based on an orthogonal sequence of polynomials, the error analysis of this measurement is 
greatly simplified; however, the rather natural interpretation of Av, Ac and as corresponding 
to structures at different densities has been lost in this scheme. 

This technique has been used in the analys is of the WiggleZ g alaxy redshift survey, the results 



of which are published in a forthcoming work (j James et al.ll2012l ). An important consideration in 



this programme is the density resolution required for a given level of accuracy in the decomposition: 
measuring higher Hermite mode coefficients requires that the genus curve be well sampled in u. 
While individual genus measurements are relatively inexpensive to compute, la rge gains are r ealise d 



when GPU computation is used, a topic explored in contemporary work bv I Wales et al.l (|2012l ). 
Analysis of this sort is well-matched to immediate use in the current generation of surveys, extending 
their capacity for description of non-Gaussian cosmic structure. 

The aspiration of this research programme is to provide an analytical account of the relation- 
ship between the physics behind non-Gaussian fields and the topological properties of the cosmic 
structures that arise as a results. Because the theoretical understanding of the relationship between 
some physical phenomena, particularly galaxy bias and nonlinear regime redshift-space distortions, 
is not established at this point, the methods presented in this paper cannot yet be exploited in 
full. Table [1] summaries this situation by describing these phenomena through the distribution 
of Hermite modes that result from the decomposition. It would be reasonable to expect that, in 
time, the impact of redshift-space distortions, galaxy bias and perhaps more general models for 
primordial non-Gaussianity, can display more complex behaviour under the decomposition, which 
would improve the prospect of understanding them through the topology of large-scale structure. 
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Fig. 7. — Same measurement as Figures [5] and [H for alternative configurations of the initial bispec- 
trum. In contrast to the local configuration, no scale-dependent bias is apparent in these results. 



Effect 


Distribution 


Comments 


Systematic selection effects 






Sampling systematic 


Delta (n = 1) 




Finite-pixel size effect 


Geometric 




Edge effects 


Uniform 




Physical phenomena 






Nonlinear gravitational evolution 




weakly nonlinear regime 


Redshift-space distortions 


Uniform 


linear regime 


Primordial non-Gaussianity 


Uniform 


assuming scale-dependent bias 


f{R) modified gravity 


Uniform 




Dark energy density 


Uniform 





Table 1: Distributional summary of results: interpreting the normalised spectrum of Hermite coeffi- 
cients On as a probability distribution, each of the physical and systematic selection effects discussed 
in this work, correspond to distributions that act multiplicatively on the Hermite spectrum. The 
current theoretical description of the effects of linear regime redshift-space distortion is that it alters 
the amplitude of the genus curve, so that the distribution of Hermite coefficients is uniform; other 
physical phenomena, with the exception of nonlinear gravitational evolution, act on the genus curve 
through the power spectrum and/or comoving volume element, though the short scale impact of 
these effects has not yet been studied. 
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Fig. 8. — Assuming a one-parameter linear relationship between scale and 02 (/nl)/«2 (gauss), 
the likelihood distribution of the gradient m is shown for each of the /nl values studied in this 
work (solid curves, main plot). Assuming that m is a normally distributed variable for each /nL; 
it is apparent that the genus statistic is not symmetrically sensitive to /nl at the 2a level for 
= {100,250}, and at the Aa level for |/nl| = 500 (inset). 
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A. Consolidation with prevailing statistics 

This Appendix presents a partial reconciliation between the formalism present ed here and 



prevailing derived statistics of genus measurements. Beginning with lPark et al.l (|l992l ). the degree 
to which a measured genus curve departs from the standard form of Equation ([2]) has been quantified 
by three numbers, termed "genus-related statistics" , that are special cases of the function 

G(a,fe,m)^ ; (Al) 

in particular 

Azy = G(-l,l,l); ylc = G(1.2, 2.2,0); = G(-2.2, -1.2, 0) (A2) 

are respectively the genus shift, cluster overabundance and void overabundance parameters. It is 
unfortunate that a little effort is necessary to recast this standard set of non-Gaussian descriptors 
in terms of the present work. A slight modification permits exploitation of a Hermite function 
recurrence relation 

oo 

^^^^l'g(y)dv = -e-^'l' ^'I'n-M, (A3) 



suggesting 

G'ia, b, m) = 4-^^ — , (A4) 



This preprint was prepared with the AAS lATpjX macros v5.2. 



-24- 



so that 



oo 




2.2 



.n=l 



V2 



i/=1.2 



(A5) 



and similarly for the void overabundance. An expression of At/' might also be derived, though it is 
regrettably difficult to reconcile the two sets of descriptors. Nevertheless, the study of deviations 
from Gaussianity across density ranges corresponding to clusters, filaments or voids is recast as 
a straightforward evaluation of Hermite functions at the density threshold values that are to be 
examined; it is stressed that the modification required to facilitate this has no bearing on the 
underlying physics. 



